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ABSTRACT 

We present an analysis of a highly asymmetric cluster merger from a Chandra observation of Abell 85. 
The merger shows significant disruption of the less massive subcluster from ram pressure effects. Nev- 
ertheless, a cold core, coincident with the cD galaxy, is observed to persist in the subcluster. We derive 
dynamical information from the motion of the cold core through the main cluster's ICM. Multiple deriva- 
tions of the velocity of the core suggest a Mach number of A4 « 1.4, or w ~ 2150 km s~^, though with 
substantial uncertainty. We construct a consistent kinematic model for the merger based on this dynam- 
ical analysis. As has been found for other such "cold fronts," conduction appears to be suppressed across 
the front. Thermal conduction may be suppressed by a magnetic field with a significant component 
perpendicular to the subcluster's direction of motion. The effect of the merger interaction in creating 
and shaping the observed radio sources is also discussed. It appears most likely that the radio source is 
due to distorted and detached lobes from the subcluster cD galaxy, rather than being a radio halo. 



Subject headings: cooling flows - 
magnetic fields 

1. INTRODUCTION 



galaxies: clusters: individual (Abell 85) — 
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Mergers of clusters of galaxies are highly energetic 
events, releasing a total kinetic energy of ^ 10^^ ergs 
into the intracluster medium (ICM). When clusters merge, 
shocks are driven into the ICM, dissipating the kinetic en- 
ergy of the merger and heating the gas. These shocks 
also have nonthermal effects, including the generation of 
turbulence in the ICM and acceleration of charged parti- 
cles to relativistic, or cosmic ray, energies. Observations 
with Chandra of merging clusters have provided new in- 
sights into the cluster merger process, including the un- 
predicted discovery of the persistence of cold cores from 
pre-merger cooling flows well into the lifetime of a merger 
("cold fronts:" Markevitch et al. 2000; Vikhlinin et al. 
2001b). 

Abell 85 is in the early stages of merging with two sub- 
clusters, each much less massive than the main cluster. 
One subcluster is merging from the southwest while the 
other subcluster is merging from the south. The south 
subcluster will be the focus of our discussion here, while 
the other subcluster and its associated radio relic will be 
discussed in a later paper. Abell 85 is unusual in being 
one of the few clusters known to be in the process of a 
merger while maintaining a moderate (107 Mq yr~^; Peres 
et al. 1998) cooling flow. Presumably, this implies that the 
merging subclusters have not yet penetrated the inner few 
hundred kiloparsecs of the cluster and have therefore not 
yet been able to disrupt the cooling flow. 

The south subcluster is more massive than the south- 
west subcluster. There has been some uncertainty in the 
past as to whether or not this subcluster is in fact merging 



with the main cluster or is merely seen against the main 
cluster in projection. Using data from ASCA, Markevitch 
et al. (1998) determined that the temperature in the re- 
gion of the subcluster is the same as or slightly greater 
than that of the rest of the main cluster at the same ra- 
dius. If the subcluster were not merging and were only seen 
it projection, its smaller mass would give it a lower tem- 
perature than that of the main cluster. Thus, the higher 
temperature indicates that the subcluster is almost cer- 
tainly interacting. 

The redshifts of the galaxies in the southern subcluster 
are slightly larger that those of the main cluster (Beers et 
al. 1991; Durret et al. 1998). This suggests that the south- 
ern subcluster is either a background cluster or that it is 
slightly in front of the main cluster and its excess redshift 
comes from its peculiar motion as it falls into the main 
cluster. Based on the analysis of Markevitch et al. (1998) 
and the observations presented in the present paper, we 
believe that it is merging with the main cluster. Thus, we 
will assume that the southern subcluster is at essentially 
the same distance as the main cluster, and that any differ- 
ence in their observed redshifts is caused by their relative 
motion along the line of sight as they merge. 

The Chandra observation and basic data reduction are 
discussed in § 2. The X-ray image is presented in § 3.1. 
In § 3.2, we analyze the spectra of interesting regions as- 
sociated with the southern subcluster. The profiles of the 
X-ray surface brightness and temperature within the sub- 
cluster and in the region ahead of the subcluster are ex- 
tracted in § 3.3. We discuss the evidence for a merger 
and X-ray determinations of the merger Mach number in 
§ 4. The pressure increase at the cold front and proper- 
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ties of the bow shock are used to derive the merger veloc- 
ity in § 4.1 and 4.2. We construct a consistent kinematic 
model for the merger in § 5. The suppression of conduction 
across the cold front is discussed briefly in § 6. There have 
been claims of a possible radio relic in this cluster as well 
(Bagchi et al. 1998), which we discuss in § 7. Our results 
are summarized in § 8. We assume Hq = 50 km s"^ Mpc^^ 
and go = 0.5 throughout this paper. At the cluster red- 
shift of z = 0.0538, 1" corresponds to 1.43 kpc. AU of the 
errors quoted are at the 90% confidence level. 

2. OBSERVATION AND DATA REDUCTION 

Abell 85 was observed with ACIS-I detector on Chan- 
dra in a single 39 ksec observation. Using the count rate 
in the S3 chip, we excluded data during two small back- 
ground flares using the lc_clean} routine written by Maxim 
Markevitch. This left 36,587 s of useful exposure time. 
Although the observation included the four ACIS-I chips 
and the S3 and S4 chips, the analysis presented here will 
be based on the ACIS-I data only. The focal plane tem- 
perature during the observation was —120 C. A raw image 
of the entire ACIS-I detector in the spectral band 0.3-10 
keV is presented in Figure 1. 
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Fig. 1. — Raw X-ray image of Abell 85 in the 0.3-10 keV band, 
uncorrected for background or exposure. All four ACIS-I chips are 
shown; the regions of reduced exposure are the interchip gaps. The 
center of the main cluster and cooling flow are located on the upper 
left chip, the center of the southwest subcluster is located just above 
the lower edge of the upper right chip, and the south subcluster is 
near the bottom and overlaps the two lower chips. 

The ACIS-I suffers from enhanced charge transfer in- 
efficiency (CTI) caused by radiation damage early in the 
mission. We corrected for the quantum efficiency non- 
uniformity and gain variations caused by this damage, but 
the degradation in spectral response was not corrected for 
due to the lack of availability of appropriate response ma- 
trices. We used version 1 of the January 29, 2001 gain file 
and response files. Because of uncertainties in the spec- 
tral response at low energy, we limit our spectral analy- 
ses to the range 0.7-10.0 keV, excluding the 1.8-2.2 kcV 
band around the mirror iridium edge. We constructed 
blank sky backgrounds from the March 23, 2001 versions 

^ see http://hea-www.harvard.edu/~maxim/axaf/acisbg/ 



of Maxim Markevitch's ACIS-I background photon lists 
using the routine make-acisbg^ . 

3. X-RAY PROPERTIES OF THE SOUTHERN SUBCLUSTER 

3.1. X-ray Image 

Figure 2 shows an adaptively smoothed image of an ap- 
proximately 6f4x6.'4 region around the southern subclus- 
ter. The image was smoothed to a signal-to-noise ratio of 
three for each smoothing beam. The same set of smooth- 
ing kernels were used to smooth the blank-sky background 
image and the exposure map. The smoothed background 
image was subtracted from the smoothed subcluster im- 
age, and the result was divided by the smoothed exposure 
map. 




Fig. 2. — Adaptively smoothed image of the south subcluster, 
corrected for background and exposure. The cold core is at the NW 
corner of the subcluster. 

The subcluster is roughly cone-shaped (Figure 2), with 
a high surface brightness knot at the northwest corner. 
This bright region is spatially extended, with a diameter 
of about 26", and is centered on the cD galaxy which is the 
brightest galaxy in the southern subcluster (see Figure 6 
below). The northern edge of this knot shows an abrupt 
surface brightness edge. This suggests that it is either a 
merger shock or a "cold front," the leading edge of a cold 
core. There is a curved tail of brighter X-ray emission 
extending to the southeast of the bright knot. 

The overall geometry indicates that the ICM in the sub- 
cluster has been affected by ram pressure from the gas 
in the main cluster. The morphology suggests that the 
south subcluster is in the early stages of merging with the 
main cluster, and that it is falling into the main cluster for 
the first time from the south. The sharp edge at the top 
of the bright knot is symmetrical about a position angle 
of —15° (15° west of north). The bulk of the subcluster 
forms a conic distribution centered about a position angle 
of 123° (58° east of south). This difference between 123° 
and 165° — (180° — 15°) may indicate that the cD galaxy 
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Table 1 
Spectral fit parameters 









Z 






Region 


model 


(keV) 






d.o.f. net counts 



cold core mekal 2.3+q;4 0.48^0.27 16^7 16 416 

cold core mckal + main cluster 2.lt°-^ O-SSlpgl 15.9 16 416 

subcluster - cold core mekal ^-^to.t 0.34^°;^^ 254.0 196 8486 

subcluster — cold core mekal + main cluster 5.5^Q g 0.36^Q'^g 254.6 196 8486 

main cluster mekal ^-Otl.t ^■^'^-olo ^^^.S 168 4527 



is moving relative to the center of mass of the subcluster, 
or that the shape has been affected by the density struc- 
ture in the subcluster and main cluster gas. In any case, 
it appears that the transverse component of the velocity 
of the subcluster relative to the main cluster lies at a posi- 
tion angle between —60° and —10°. The position angle of 
the center of the main cluster from the subcluster is about 
+ 13°. This implies that this is an offset merger; the col- 
lision is occurring with a nonzero impact parameter and 
angular momentum. 

Several point X-ray sources are also seen in the region 
of the subcluster (Figure 2). Only two of these have a an 
optical counter part: the source at R.A. = 00'^41™50!4, 
Dec. = — 9°25'48" is coincident with the nucleus of the 
galaxy PGC 93226, which is a cluster member and a ra- 
dio source (radio source B in Figure 6 below); and the 
source at R.A. = 00*^41™59!0, Dec. = -9°24'49" is coinci- 
dent with a very faint galaxy of unknown redshift (source 
1-2117; Slezak et al. 1998). 

3.2. X-ray Spectra 

We extracted the X-ray spectra of the bright knot co- 
incident with the cD galaxy at the top of the subcluster 
("cold core"), of the remainder of the subcluster, and of 
the main cluster gas at a similar projected distance from 
the center of the main cluster. The spectrum of the bright 
knot was extracted from the elliptical region at the center 
in Figure 3a. The spectrum from the remainder of the 
subcluster was extracted from a polygonal region encom- 
passing most of the rest of the subcluster, minus the point 
sources within that region. The main cluster spectrum 
came from an annular pie wedge to the east of the sub- 
cluster and which encompassed the same radii from the 
cluster center as the subcluster. The spectral fits in these 
regions are presented in Table 1. The spectra were grouped 
to have a minimum of 20 counts per channel. As described 
in § 2, we restricted our spectral analysis to the range 0.7- 
10.0 keV, minus the band from 1.8-2.2 keV. In the fit to 
the spectrum of the main cluster, we cut the spectrum oS 
at 9.0 keV because the spectrum of this diffuse emission is 
dominated by background above this energy. The spectra 
were fit within xspec using the mekal model for the thermal 
emission. The absorption column was fixed at the Galactic 
value of 2.85 x lO^" cm'^ (Dickey & Lockman 1990). Since 
the observed emission from the cold core and subcluster 
presumably contain emission from the main cluster seen 
in projection in front of and behind these regions, we also 
fit the spectrum of the cold core and subcluster with two 
mekal thermal components, with the shape of the hotter 
component fixed at the values found for the main cluster, 



and the normalization determined by the relative areas of 

the regions. 

We find that the temperature of the bright knot at the 
top of the subcluster (2.1^0^ keV) is much lower than the 
temperature of the remainder of the subcluster or of the 
surrounding gas from the main cluster. This shows that 
this X-ray bright and dense region is a cooling core associ- 
ated with the central region and central cD galaxy in the 
subcluster. The sharp surface brightness discontinuity at 
the northern edge of this knot must be a cold front, rather 
than a merger shock, since the compressed gas has a lower 
temperature and specific entropy than the less dense gas 
(Markevitch et al. 2000; Vikhlinin et al. 2001b). 

We also fit the cold core with a cooling flow model, with 
and without an additional foreground and background 
contribution from the main cluster. In both cases, the 
gas was allowed to cool to the minimum allowable temper- 
ature, essentially zero. For the fit without the additional 
component for the main cluster, we fixed the maximum 
temperature and abundance to those from the third fit in 
Tablet. We found a cooling rate of 7.3t2;^ Mq yr"!. With 
an additional mekal model component set to the main clus- 
ter parameters, and with the maximum temperature and 
abundance set to the parameters from fit number 4 in Ta- 
ble 1, we derived a cooling rate of 6.8±0.6 Mq yr^-'^. Both 
fits are consistent with a low present cooling rate. Both fits 
also had a significantly worse reduced than did either 
of the first two models presented in Table 1. 

We fit a single temperature model to the spectrum of the 
subcluster minus the cold core and found a temperature of 
6.31q 5 keV. If we add a model component for the emission 
from the main cluster, we find a subcluster temperature 
of 5.5^Q g kcV. The former temperature is slightly lower 
than the value given by Markevitch et al. (1998), but is 
consistent to within the errors. Our extraction region is 
smaller than that used by Markevitch et al. (1998) due to 
the much poorer angular resolution of ASCA. Thus, the 
ASCA spectrum may have included more emission from 
the main cluster, which is hotter. In any case, the bulk 
of the subcluster is hotter than might be expected for a 
cluster of this mass and X-ray luminosity, which may indi- 
cate that the much of the subcluster gas has been heated 
by shocks or adiabatic compression associated with the 
merger. 

3.3. Temperature and X-ray Surface Brightness Profiles 

We measured the temperature and surface brightness 
gradients inside the subcluster and in front of the cold front 
(Figure 3). The temperature measurements in front of the 
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Fig. 3. — (a) Elliptical annular wedge regions used for determining temperatures are shown superposed on a gaussian smoothed image 
of the south subcluster. (b) Temperature profile in the regions shown in (a). The dotted line indicates the position of the cold front, with 
positive radii in the direction of motion of the subcluster (i.e., ahead of the cold front). The dashed lines indicate the temperatures from 
the single-temperature fits given in Table 1. The fit to the subcluster is shown for the radii over which it was determined, (c) X-ray surface 
brightness from 0.3-10 keV in a set of elliptical annular wedges with the same shape and orientation as those in (a), but more closely spaced. 
(d) X-ray surface brightness values from (c), after correction for the foreground and background emission from the main cluster. The dashed 
curve in ('cj is the projected surface brightness contribution from the main cluster, which is used to correct the values in (d). 



cold front were made by extracting spectra in a wedge 
of elliptical annuli whose curvature matched that of the 
cold front. The measurements within the subcluster were 
made from spectra also accumulated from annular wedges 
using ellipses self-similar to those in front of the cold front. 
The size and orientation of the wedges inside the subclus- 
ter were determined by the edges of the subcluster, and 
therefore were not oriented 180° from the wedge in front 
of the subcluster. While the annular wedges in front of the 
subcluster were centered along a line 15° west of north, the 
wedges inside the subcluster were centered on a line ^ 58° 
east of south. The regions used to extract the spectra are 
shown in Figure 3a. 

We fit single-temperature models to these spectra, with 
the absorption column set to the Galactic value. The re- 
sulting temperatures from fits to these spectra are shown 
in Figure 3b. The gas in front of the cool core is hot, and 



the temperatures are consistent (within the large errors) 
with the temperature in the main cluster at this radius. 
The temperature in the cool core of the subcluster is quite 
low. Behind the cool core, the temperatures in the sub- 
cluster rise up to moderately high values. 

X-ray surface brightness measurements were made in el- 
liptical annular wedges with the same shape as those used 
to extract the spectra, but with smaller widths. The re- 
sulting surface brightness profile in the 0.3-10 keV band is 
shown in Figure 3c. As is clear from the image (Figure 2), 
the highest surface brightness is associated with the cool 
core. There is a very sharp surface brightness discontinu- 
ity (a factor > 5) at the northern edge of the cool core. 
The combination of the surface brightness discontinuity 
with the low temperature in the bright region shows that 
this is a cold front. The subcluster south of the cool core 
is much brighter than the gas ahead of the cold front. 
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If one assumes that the subcluster is merging for the 
first time, that the motion is transonie, and that the mass 
of the subcluster is much smaller than that of the main 
cluster, one would expect the gas far ahead of the cold 
front would be undisturbed main cluster ICM. Also, main 
cluster emission may be projected in the foreground and 
background of the subcluster. To determine the contri- 
bution of undisturbed main cluster emission, we measured 
the surface brightness profile of the main cluster in a wedge 
to the southeast; this region is essentially the same region 
as that occupied by the subcluster, but refiected across 
the north-south axis of symmetry of the main cluster. We 
fit a /3-model to the surface brightness at projected radii 
from 380 to 680 arcsec, which covers the range of radii 
containing the subcluster. 

The main cluster surface brightness determined from 
this fit is shown as a dashed curve in Figure 3c. More 
than ~40 kpc ahead of the cold front, the X-ray surface 
brightness is consistent with the undisturbed main cluster 
emission within the errors. However, there is some evi- 
dence for a rise in the surface brightness just ahead of the 
cold front; the four values at ~0-40 kpc are all slightly 
higher than expected. This may indicate that the main 
cluster gas is compressed ahead of the subcluster and cold 
front, by a bow shock and/or by adiabatic compression. 

To show more clearly the excess X-ray emission asso- 
ciated with the subcluster and any compression of main 
cluster gas ahead of the cold front, in Figure 3d we sub- 
tract the fit to the undisturbed main cluster emission from 
the surface brightness values in Figure 3c. In most of the 
region in front of the subcluster, there are only upper lim- 
its on the excess emission. There may be some excess 
emission just ahead of the cold front, but in the residual 
surface brightness profile only the point from ~1-10 kpc 
appears to be significantly increased. In the residual pro- 
file, the surface brightness of the subcluster is relatively 
uniform except for the brighter cool core. Thus, the ap- 
parent fall-off in the surface brightness of the subcluster 
with increasing distance from the cold front in Figure 3c 
may actually be due to projected main cluster emission. 

We determined the gas densities in the regions aroimd 
the cold front and subcluster by deprojection. We as- 
sumed diS'erent geometries for the subcluster and for the 
gas ahead of the cold front. For the densities inside the 
subcluster, we assumed the subcluster to be a cone open- 
ing up behind the cold front, with an opening angle de- 
termined from the image. We further assumed that this 
cone's axis of symmetry lies in the plane of the sky. We de- 
termined the gas density well ahead of the cold front from 
the /3-modcl fit to the surface brightness of the main cluster 
discussed above. We also assumed a spherically symmetric 
main cluster to do the deprojection. We used the depro- 
jected gas densities and the temperatures from spectral 
fits to determine the pressures across the cold front. 

4. HYDRODYNAMICAL ANALYSIS OF MERGER 

We now use the gas temperatures, densities, and pres- 
sures derived in § 3.3 to analyze the kinematics and hydro- 
dynamics of the subcluster merger. Our treatment closely 
follows that in Vikhlinin et al. (2001b). A schematic view 
of the geometry of the fiow of main cluster gas near the 
cold front is shown in Figure 4), which is adapted from 



Vikhlinin et al. (2001b). A bow shock will be present 
ahead of the cold front if the merger velocity of the cold 
front relative to undisturbed main cluster gas is super- 
sonic. (There is also a shock in the cold front, but the 
shocked region will be very narrow if the cool core is much 
denser than the main cluster gas.) If the motion is sub- 
sonic, the flow around the cold front is continuous. For 
any blunt cold front, there will be a point in front of the 
cold front where the velocity of the main cluster gas is 
zero. This stagnation point is labeled "st" in Figure 4). 

4.1. Stagnation Pressure at Cold Front 

The ratio of the pressures in the main cluster gas far 
ahead of the cold front to that at the stagnation point, 
combined with temperature measurements of the gas, 
can be used to determine the velocity of the cold front 
(Vikhlinin et al. 2001b). Ideally, we would measure the 
pressure at the stagnation point in the hot gas. However, 
since the surface brightness of the hot gas at the stagna- 
tion point is actually quite low, we instead measure the 
pressure in the cool core just behind the stagnation point 
where the surface brightness is much higher. Because the 
cold front is a contact discontinuity and the gas is moving 
subsonically near the stagnation point, the pressure across 
the cold front is expected to be continuous. The pressure 
difference between the stagnation point and a point far 
upstream (region 1 in Figure 4) must be caused by com- 
pression of the gas in region 2 by a bow shock and/or 
adiabatic compression. 




bowshock? 



Fig. 4. — Schematic diagram of the flow of hot main cluster gas 
around a blunt cold front. A bow shock will be present ahead of 
the cold front if the merger is supersonic. The stagnation point is 
labeled "st" . Region 1 is unshocked gas of the main cluster; region 2 
is gas which has passed through the bow shock (if present). 

The ratio of the pressure at the stagnation point to the 
pressure in the far upstream region 1 is given by (e.g. Lan- 
dau & Lifshitz 1959, §114) 



Pst^ 

Pi 



(1 + 21^^2).-! 



M<1, 
M>1. 



(1) 
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Here Pst and Pi are the pressures at the stagnation point 
and in region 1, respectively and 7 = 5/3 is the adiabatic 
index for a fully ionized plasma. A4 = Vi /Cgi is the Mach 
number of the cold core, vi is the cold core's velocity rel- 
ative to the upstream gas, and Cgi is the sound speed in 
that gas. 

We do indeed measure a higher pressure inside the cold 
core as compared to the undisturbed gas in front of the 
core. The best fit measurement is Pst /Pi = 3.4 which im- 
plies a Mach number of 1.4. This mcasiircment assiimes 
that the separation between the main cluster and the sub- 
cluster is equal to their projected separation, so Pi is mea- 
sured at the projected radius of the subcluster within the 
main cluster. The formal errors yield a wide range in al- 
lowable Mach numbers, from to 3.3, but since not all the 
sources of error are independent, the actual error is prob- 
ably somewhat smaller. While a Mach number of zero 
is allowed, the morphology of the system makes such a 
value highly unlikely. In any case, the best fit value of the 
pressure ratio requires that the merger motions be slightly 
supersonic. 



4.2. Possible Bow Shock 

For a supersonic cold core, a bow shock should also form 
in front of the cold core. Assuming that the density and 
temperature of the gas in region 1 is constant, the bow 
shock should have a predictable "stand-ofF" distance, ds, 
which is the shortest distance from the stagnation point 
to the bow shock (Vikhlinin et al. 2001b). This distance 
can be calculated using the approximate method given by 
Moekcl (1949), and depends only on the value of A4 and 
on the shape of the cold front. A useful plot of the stand- 
off distance versus Mach number is given in Sarazin (2002, 
Figure 4). For > 2, the stand-off distance of the shock 
is not very sensitive to the value of M, while for smaller 
Mach numbers the distance increases rapidly. For the best 
fit value of = 1.4, the expected stand-off distance is 
rfg ~ 18 kpc if we treat the cold front as a spherical sur- 
face with a radius of curvature of 19 kpc. Such a bow shock 
would compress the gas in region 2, and should produce 
a measurable increase in the X-ray surface brightness, Ix , 
in that region. A possible surface brightness excess is seen 
in the ~20 kpc immediately upstream from the cold front, 
but it is significant at only slightly greater than the 1.7a 
level. While this surface brightness excess is not visible in 
the image of the cluster, we might expect such a feature to 
be more visible in profile given the significant azimuthal 
averaging done to create the profile. 

Given enough source photons, the spatial resolution 
of Chandra would be sufficient to measure the expected 
stand-off distance of the bow shock, which corresponds to 
~5-15". However, the surface brightness in the hot gas 
ahead of the cold front is too low to allow dg to be accu- 
rately determined from the available data. Thus, all we 
can conclude is that the expected values of the stand-off 
distance are consistent with the (marginal) evidence for an 
increase in the X-ray surface brightness within ~20 kpc of 
the cold front. 

We can also use the Rankine-Hugoniot shock jump con- 
ditions at the putative bow shock to independently de- 
termine the Mach number (e.g. Landau & Lifshitz 1959, 



§85). The shock jump conditions yield 



C 7 + IAI2 ^ ^+1' V ^ 

where C = P2/P1 is the shock compression. Because we 
do not have spectra or temperatures determined on the 
scale of the bow shock, we estimate the shock compres- 
sion from the small increase in the surface brightness as 
C « [1x2/ IxiY^'^ ■ The observed surface brightness in- 
crease in the first 10 kpc is a factor of ^ 1.7, which im- 
plies that C ~ 1.3. This imphes M ^ 1.2. It is likely 
that the finite resolution with which the surface bright- 
ness contrast was determined and projection effects cause 
the shock compression to be underestimated. Thus, this 
value for the Mach number is consistent with that deter- 
mined from the pressure increase at the stagnation point. 
If we assume the Mach number determined by the stagna- 
tion condition M. ^ lA, the expected shock compression is 
C ~ 1.6. Projection effects (if the cluster is not moving in 
the plane of the sky) would cause us to overestimate the 
stand-off distance and underestimate the Mach number. 
Projection could also cause us to inaccurately determine 
the true shape of the cold front. 

5. MERGER KINEMATICS 

Since it provides a consistent fit to the stagnation pres- 
sure, the bow shock compression, and the bow shock 
stand-off distance, we will adopt the merger Mach num- 
ber oi ^A K. 1.4. The sound speed in the upstream gas is 
» 1540 km s~^, so a Mach number of 1.4 implies a merger 
velocity of u « 2150 km s~^. 

5.1. Kinematic Model 

We now construct a kinematic model for the merger 
which is consistent with the X-ray and optical observa- 
tions of the main cluster and subcluster. Because the er- 
rors on our determination of the Mach number are large, 
the parameters of this model are not well constrained. We 
therefore do not suggest that this model is the only pos- 
sible interpretation of the data, but that it is merely a 
"toy model" that is adequate to explain the data given 
the best- fit values from the various hydrodynamic tests. 
The model, then, is presented as illustrative rather than 
interpretive, using the best-fit parameters from the hy- 
drodynamic analyses as the primary constraints of the 
model. We also discuss the implications of this model for 
the merger. 

The parameters of the model are shown in Figure 5. We 
will approximate the subcluster as a point mass with a 
single velocity relative to the main cluster. We will also 
assume that the mass of the subcluster is small relative to 
that of the main cluster, so we can treat the subcluster as 
a test particle falling into the extended mass distribution 
of the main cluster. We put the center of the main clus- 
ter at the center of our coordinate grid, and using polar 
coordinates, we define the a;-axis to be parallel to North 
and the z-axis to be the line of sight, with the positive 
^;-axis extending away from the observer. Let the vector d 
be the position the subcluster relative to the main cluster 
(direction from the main cluster to the subcluster). The 

components of d are defined by its magnitude d, the angle 
to the line of sight 64 (from the positive 2;-axis), and the 
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position angle on the plane of the sky (pd measured coun- 
terclockwise from the north. Similarly, let v be the velocity 
of the subcluster relative to the main cluster, with magni- 
tude V, and direction given by the angles 9y and (j)y. We 
also define tjj to be 180° minus the angle between d and v. 

On the X-ray image (Figure 1), the subcluster center 
is located at a position angle of <pd = 194° from the main 
cluster. The position angle of the direction of motion, (jj^ is 
less certain, since the curvature of the cold front suggests 
~ —15° whereas the body of the subcluster is more 
consistent with ^ —58°. We will adopt the average 
value of « —36°. 

We can estimate the radial component of the relative 
velocity of the subcluster from the optical rcdshifts of the 
main cluster and subcluster. We will adopt the velocity of 
the subcluster cD galaxy as representative of the subclus- 
ter; in any case, it is most closely related to the cold front, 
which was used to derive the merger velocity. Beers et 
al. (1991) give the redshift of the subcluster cD galaxy as 
z = 0.05633±0.00012, while the galaxies within 5'.25 of the 
center of the cluster have a mean redshift of 0.0538±0.0050 
(Durret et al. 1998). Combining the line-of-sight veloc- 
ity determined from the optical redshifts with the velocity 
merger velocity determined in § 4 from the X-ray data, we 
find that the direction of motion of the subcluster relative 
to the main cluster is between 9-^ = 63° and 76° from the 
line of sight, or 14° to 27° from the plane of the sky. We 
adopt the average value of ^„ « 71°. 

a? 




Fig. 5. — Schematic diagram of the kinematic model. The position 
of the subcluster relative to the main cluster is given by the vector 
d, and v is the velocity of the subcluster relative to the center of the 
main cluster. The z-aKis is along the line-of-sight. The dashed lines 
arc projections of the direction and velocity vectors onto the plane 
of the sky. The dotted line is parallel to the N-axis. The small inset 
shows the projection of the direction and velocity vectors onto the 
N-z plane, with the dotted line parallel to the 2-axis. 

It is difficult to determine the angle between the separa- 
tion of the main cluster and subcluster, d, and the line-of- 
sight. Initially, we will assume that the main cluster and 
subcluster are separated in the plane of the sky {64 = 90°), 
so that the separation between the main cluster and clus- 



ter is equal to the projected separation, d = 730 kpc. The 
angle between d and v is then = 53°, and the impact 
parameter of the merger is 6 w 580 kpc. 

5.2. Infall Velocity 

We now compare the merger velocity inferred from the 
hydrodynamics of the cold front and subcluster with the 
infall velocity expected for the subcluster and main clus- 
ter. Reiprich & Boringer (2002) derive a virial mass for 
the main cluster of M200 = 1.080 x 10^^ Mq for a virial 
radius r2oo = 2.66 Mpc. We will assume that the main 
cluster mass is much larger than that of the south sub- 
cluster. We will assume that the subcluster has fallen into 
the main cluster from its turn-around distance of 5.5 Mpc, 
which is the value if the age of the Universe is 13 Gyr. As 
noted above, we initially assume that the main cluster and 
subcluster are separated in the plane of the sky. 

We will consider two models for the mass distribution 
and potential of the main cluster. We first consider a sin- 
gular isothermal sphere out to the virial radius, which has 
a potential given by 

$(r) = -2a2x( 1 + MWr), r ^ r^oo , 
^ ^ \ r20o/r , r > r2oo • ^ ^ 

Here, a = (GM2oo/2r2oo)^/^ « 934 km s"^ is the velocity 
dispersion. The infall velocity at the subcluster's current 
projected distance of 730 kpc is 2520 km s~^. We also 
consider a model in which the density within the virial ra- 
dius is given by the Navarro, Frenk, & White model (1997, 
hereafter NFW), for which the potential is 

[ J , r> r2oo • 

Here, is the scale radius, and x = r/rg. We adopt a 
concentration parameter c = r^oo/vs = 10, which is con- 
sistent with NFW's simulations for cluster-mass halos. For 
this potential, the predicted infall velocity at the projected 
distance is 2740 km s~^. 

For either potential, the infall velocity at the projected 
separation is somewhat larger than the velocity we deter- 
mined from the X-ray observations of the merger hydrody- 
namics. Given the large errors in the determinations of the 
velocities, this difference may not be significant. As first 
noted by Markevitch et al. (1999), the degree of agreement 
between the merger velocity determined by hydrodynam- 
ical measurements and that expected from infall can be 
used to test the hypothesis that the intracluster medium is 
predominantly a non-relativistic, thermal plasma. That is, 
the calculation of of the merger velocity from shock condi- 
tions (equation 2) assumes that the merger shock energy is 
thermalizcd, and is not converted into relativistic particles, 
or magnetic fields, or turbulence. All of the hydrodynamic 
diagnostics require that the intracluster medium act as a 
7 = 5/3 gas. Thus, the difference between our hydrody- 
namical estimate of the merger velocity and the predicted 
infall velocity (assuming the main cluster and subcluster 
are separated in the plane of the sky) could suggest that 
the kinetic energy of the merger is not thermalizcd par- 
ticularly efficiently, but instead goes partially into turbu- 
lence, magnetic fields, or relativistic particles. As noted by 
Markevitch et al. (1999), this argument is somewhat cir- 
cular, as the mass of the cluster was also determined from 
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hydrostatic equilibrium assuming purely thermal pressure 
support. Given the uncertainties in the determination of 
merger velocity and infall velocity, we will instead take the 
crude agreement between the two speeds as an indication 
that at least a significant fraction of the merger energy 
(> 50%) is thermalized. 

Our initial estimate of the infall velocity was based on 
the assumption that the main cluster and the subcluster 
were separated in the plane of the sky (6*^ = 90°). If this 
is not true, the actual separation d will be larger than 
the projected separation of 730 kpc, and the predicted in- 
fall velocity will be lowered. To illustrate this effect, we 
construct a consistent model for the merger geometry and 
kinematics in which the merger velocity equals the pre- 
dicted infall velocity. As wo move the subcluster further 
out in the main cluster potential, the density and hence the 
pressure in the ambient medium drops, thereby increasing 
the pressure ratio used to determine the Mach number. 
This in turn increases the Mach number we would measure, 
lessening the need to place the subcluster significantly in 
front of or behind the main cluster. 

Given the large errors in the two numbers, our model 
is certainly not a unique solution, but is consistent with 
the current best-fit values of the parameters. For the 
isothermal potential (equation 3), this consistent solution 
requires that d « 820 kpc, while for the NFW potential 
(equation 4) we find d w 860 kpc. The Mach numbers 
we derive arc A4 = 1.6 for the isothermal potential and 
M. = 1.7 for the NFW potential. The corresponding ve- 
locities are v « 2460 km s~^ and v « 2610 km s~^ re- 
spectively. For both solutions, we assume that the sound 
speed does not vary over the range of radii in question. 
If we were to allow the sound speed to vary, it would de- 
crease slightly at larger radii, increasing the physical sep- 
aration that we determine. To be specific, we will adopt 
the isothermal result. This distance implies that the an- 
gle between the separation and the line of slight is either 
9d w 64° or 9d ~ 116°. The former value implies that 
the ~ 66°, which means that the subcluster is moving 
nearly perpendicular to the radius from the center of the 
main cluster, and about to start exiting the cluster. The 
observed morphology of the X-ray image of the subclus- 
ter seems inconsistent with this interpretation. Thus, we 
adopt the solution with Od « 144° and tp ~ 46°, in which 
the subcluster is moving into the main cluster, probably for 
the first time. In summary, our consistent kinematic model 
has V « 2460 km s-\ 0y ^ 71°, <py -36°, d « 820 kpc, 
6»d » 116°, (f)d » 194°, and V » 46°. The subcluster is 
closer to us than the main cluster, and is falling into the 
main cluster. 

5.3. Angular Momentum and Impact Parameter 

As noted above in § 3.1, the direction of the merger ve- 
locity is not parallel to the separation of the centers of the 
main cluster and subcluster. The transverse component of 
the merger velocity lies at an angle of ~50° with respect 
to the projected separation of the two clusters. This im- 
plies that this is an offset merger; the collision is occurring 
with a nonzero impact parameter and angular momentum. 
If we adopt the consistent model for the merger kinemat- 
ics which we have just discussed, the angle between the 
velocity and the separation is V ~ 46°, and the impact 



parameter for the collision is 6 = dsini/; f« 750 kpc. This 
is 3.7 times the core radius we determine from a /3-model 
fit to the cluster. 

A useful dimensionless form for the angular momentum 
is given by the A parameter, defined as (Peebles 1969) 



Here J is the total angular momentum of the merged clus- 
ter, E is its total energy, and M is its mass. We estimated 
the value of A implied by the merger velocity and impact 
parameter of the subcluster by differentiating equation (5), 
assuming the mass of the subcluster is much smaller than 
that of the main cluster. We assumed that the orbital 
angular momentum of the subcluster was parallel to the 
initial angular momentum of the main cluster, and ignored 
the initial internal angular momentum of the subcluster. 
Using the kinematic parameters for our consistent model 
for the merger, we find A » 0.21. Fixing v, 6^, and the pro- 
jected angles (/>d and at the values from the consistent 
model, and allowing the value of 9d to vary, we find a mini- 
mum value of A of about 0.16. If the initial spin of the main 
cluster is not aligned with the orbital angular momentum 
of the subcluster, the value of A would be an upper limit. 
Those values are both somewhat larger than the median 
values of 0.05-0.1 expected from tidal effects in large scale 
structure. This may reflect the large uncertainties in the 
kinematic parameters. If this large angular momentum is 
correct, it might be the result of tidal effects associated 
with the triple merger occurring in Abell 85 (i.e., the fact 
that there are two merging subclusters). Alternatively, it 
may be that mergers with small subcluster have a larger 
range of values of A, which average out when many small 
subclusters merge to form a larger halo. 

6. SUPPRESSION OF CONDUCTION ACROSS THE COLD 

FRONT 

As discussed in § 3.1, the cold front is seen as a sharp 
surface brightness discontinuity, with a dramatic increase 
in surface brightness in the cold gas over only a few kilopar- 
secs. From the surface brightness profile in Figure 3, the 
change in surface brightness occurs over at most 20 kpc — 
half the width of the elliptical region on the cold core or 
the width of about 2 bins ahead of the cold front. The 
raw image suggests that it is actually quite a bit narrower 
than this, but we are limited by photon statistics to the 
aforementioned resolution. Unfortunately, the same pho- 
ton statistics prevent us from measuring the temperature 
gradient on this size scale. However, because the pressure 
is continuous on the largo scale of our measurements and is 
presumed to be continuousl on the smaller scale of the cold 
front as wel, the observed density gradient should be ac- 
companied by a temperature gradient of the opposite sign 
and with the same length scale. Therefore, while the width 
of the gradient is determined from the surface brightness, 
we can assume that the same width applies to the tem- 
perature gradient. Thermal conduction should smear out 
any such sharp edges to a length a few times the electron 
mean free path in a relatively short time. However, con- 
duction appears to be suppressed in the case of Abell 85, 
as it docs in other clusters with observed cold fronts (e.g. 
Abell 2142, Abell 3667; Ettori & Fabian 2000; VikhUnin 
et al. 2001b). 
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Ettori & Fabian (2000) showed that thermal conduction 
should smear out the temperature gradient to a width Sr 
on a characteristic timescale 

Sv 

Sr = ^, (6) 

V 

where 

_ 2 K djkBTe) 

3 rieksTe dr 
is the characteristic velocity of the diffusion. Here, 

K = 8.2x lO^o ( Jl^) ergs-i cm"! keV-^ (8) 
\ 10 keV J 

is the thermal conductivity, and Ue and Tg are the elec- 
tron number density and temperature, respectively. For 
the upper limit on the width of the cold front of 20 kpc, 
the diffusion timescale is 2.0 x 10^ yr. 

At the subcluster's current velocity and distance from 
the cluster center, the relevant timescale for interaction 
is roughly d/v k, 5.6 x 10^ yr. This means that in order 
for thermal conduction to have failed to erase the sharp 
edge of the cold front, conduction must be suppressed by 
at least a factor of 280-2700. In fact, since the rate of 
conduction is independent of density, the time over which 
conduction has had a chance to act is probably somewhat 
longer, meaning that the degree of suppression is probably 
even higher. 

One mechanism that has been suggested for this sup- 
pression is the existence of a magnetic field perpendicular 
to the direction of diffusion, i.e. parallel to the surface of 
the discontinuity (Vikhlinin et al. 2001b). A tangled mag- 
netic field would serve the same purpose, with a maximum 
loop size equal to the width of the front. 

7. subcluster/radio source interaction 

The possible existence of a large scale radio halo or relic 
in the southern subcluster in Abell 85 was raised by Bagchi 
et al. (1998). As shown in Figure 6, the diffuse emis- 
sion was later resolved by Giovannini & Feretti (2000) into 
what may be a tailed source associated with a dead or dy- 
ing AGN in the dominant galaxy of the subcluster, where 
the nuclear emission is very faint compared to the radio 
lobes. If this is in fact its origin, the shape of the tail is 
well explained by the merger interaction. Figure 6 shows 
that the source associated with the cD galaxy in the sub- 
cluster (Source D) has possible weak nuclear emission, and 
a lumpy C-shaped extended source to the southeast of the 
cD galaxy. The radio emission has no significant extent 
north of the subcluster's leading edge, and the bulk of the 
radio emission follows the bright X-ray arc through the 
subcluster. If it is indeed a tailed galaxy, this shape im- 
plies that the same ram pressure forces which have shaped 
the subcluster have also bent the lobes of the AGN into 
their current shape. 

An alternative explanation for the diffuse radio emission 
is that this is a small, merger-induced halo in the subclus- 
ter. The fact that its surface brightness correlates roughly 
with that of the X-ray gas is consistent with findings for 
other halos, as is its steep spectral index (oq | ~ 2—2.5; 
Giovannini & Feretti 2000). On the other hand, the ra- 
dio source is quite overluminous compared to the expected 
radio power derived from the empirical relation between 



X-ray temperature or luminosity and radio halo power 
(Liang et al. 2000; Feretti 2000). For the observed X- 
ray luminosity of the subcluster of 5 x 10^'^ erg s~^ (rest 
frame 0.1-2.4 keV), the empirical relation would require 
a monochromatic radio power of only 4 x 10^^ W Hz"^ 
at 1400 MHz in the cluster rest frame, compared to the 
observed value of 1.8 x lO^^ W Hz"! (Feretti 2001) Fur- 
thermore, halos have previously only been observed in the 
hottest and most massive clusters, never in a cluster as 
small and cool as this subcluster. While previous observa- 
tions may have been biased towards finding halos in hot 
clusters (Kempner & Sarazin 2001), it would nonetheless 
be surprising to find a halo in such a cool, low luminos- 
ity cluster. Our observations are therefore more consistent 
with the interpretation of the source being a tailed radio 
galaxy (Giovannini & Feretti 2000) than with it being a 
radio halo or relic. 

Source "B" in Figure 6 is a narrow angle tail (NAT) 
source studied in detail by O'Dea & Owen (1985). It has 
a redshift of z = 0.0579 (Wegner et al. 1999), which is 
within the dispersion of the main cluster. Since the line 
of sight component of the subcluster's infall velocity is so 
small, however, it is not possible to determine whether or 
not this galaxy is a member of the subcluster or of the 
main cluster based solely on radial velocity information. 
The direction of the tail, plus the fact that the radio tail 
does not appear to be interacting with the flow around the 
subcluster, suggest that it is seen in projection in front of 
or behind the subcluster, and that it is a member of the 
main cluster. 




Fig. 6. — Radio contour map at 90 cm of the "B" and "D" sources 
(from Giovannini & Feretti 2000) in the region of the southern sub- 
cluster, overlaid on the adaptively-smoothed Chandra image of the 
subcluster (color, from Figure 2). The D source appears to be as- 
sociated with the cD galaxy in the south subcluster, while the B 
source is associated with another cluster galaxy. 



8. SUMMARY 

Our analysis of the south subcluster in Abell 85 from 
~ 37 ksec of Chandra data has revealed several interesting 
features. The most obvious is a confirmation that the sub- 
cluster is indeed merging with the the main cluster. The 
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subcluster contains a remnant cold core which has survived 
the early stages of the merger. It is smaller and more dis- 
crete than similar structures found in other clusters such as 
Abell 2142 (Markevitch et al. 2000), Abell 3667 (Vikhlinin, 
Markevitch, & Murray 2001b), and RX J1720. 1+2638 
(Mazzotta et al. 2001), and is perhaps more akin to the 
"bullet" in 1E0657-56 (Markevitch et al. 2002) or the 
"tongue" seen in Abell 133 (Fujita et al. 2002). 

Based on the ratio of the pressure at the stagnation 
point of the cold front to that far upstream, on the stand- 
off distance of a possible bow shock, and on the shock 
compression from the bow shock, we find a consistent 
Mach number and velocity for the merger oi A4 ~ lA and 
V f» 2150 km s~^. By comparing this velocity to the radial 
velocity of the subcluster relative to that of the main clus- 
ter, we have determined that the merger velocity is about 
19° from the plane of the sky. We find a consistent kine- 
matic model for the merger in which the subcluster is in 
front of and falling into the main cluster. This model is 
consistent with the expected merger velocity if the subclus- 
ter and main cluster have fallen towards one another due 
to gravity from their turn-around distance in the Hubble 
flow. 

The X-ray observations indicate that this is an offset 
merger with a finite impact parameter and a significant 
angular momentum. A crude estimate based on our con- 
sistent kinematic model suggests an angular momentum 



parameter of A ~ 0.2, which is somewhat larger than the 
median values expected due to tidal torques. 

Magnetic fields in the cold core may be responsible for 
suppressing thermal conduction across the cold front. This 
would explain the sharpness of the front, which should be 
smeared out by conduction in the absence of a magnetic 
field. The magnetic fields in the cold core may be high as 
a result of a cooling flow or the AGN located in the central 
cD galaxy. 

We confirm the assertion that the diffuse radio structure 
in the subcluster is not a cluster radio halo or relic, but is 
more likely to be a tailed galaxy with a weak or dead nu- 
cleus. We also show that its morphology has been shaped 
by ram pressure in the merger interaction. 
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